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Abstract 



In a previous work [J. Math. Phys. 35 (1994), 2539-2551], generahzed hypergeo- 
r~| \ metric functions have been used to a give a rigorous derivation of the large s asymp- 

0-[ totic form of the general /3 > gap probability £'^'^'''^(0; (0, s); /3a/2), provided both 

/3a/2 € Z>0 and 2/13 £ Z"*". It shown how the details of this method can be extended 
to remove the requirement that 2//3 S Z"*". Furthermore, a large deviation formula 
for the gap probability Ep{n] (0, x); ME_a^7v(A"^''^e^^'^/^)) is deduced by writing it in 
terms of the charateristic function of a certain linear statistic. By scaling x = s/{AN)'^ 
and taking N ^ oo, this is shown to reproduce a recent conjectured formula for 
^^^■•'^(n;(0,s);/3a/2), pa/2 G Z>o, and moreover to give a prediction without the 
latter restriction. This extended formula, which for the constant term involves the 
Barnes double gamma function, is shown to satisfy an asymptotic functional equa- 
tion relating the gap probability with parameters (/3,n, a), to a gap probability with 
parameters (4//3, n', a'), where n' = P{n + l)/2 - 1, a' = /3(a - 2)/2 + 2. 

1 Introduction 

The topic of eigenvalue spacing distributions within random matrix theory is of interest 
both for its practical utility in comparisons with experimental data, and for its rich math- 
ematical content leading to many explicit functional forms. Moreover, the function theory 
associated with the latter is intimately related to integrable systems theory, and this in 
turn offers a number of powerful asymptotic methods to further quantify the spacing dis- 
tributions. We refer to [I5l [1] for survey articles relating to explicit functional forms, and 
to [23] for a recent review on the associated asymptotics. 

Spacing distributions depend on a symmetry parameter /3, and the region of the eigen- 
value spectrum being scaled. In relation to /3, according to Dyson's three fold way |10j . 
random matrix ensembles corresponding to quantum systems with a time reversal symme- 
try T such that = 1 must have a probability density function invariant under conjugation 
by orthogonal matrices. If the quantum system does not have a time reversal symmetry. 
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or if the time reversal symmetry is such that = —1, the probabihty density function 
must be invariant under conjugation by unitary matrices, or by unitary symplectic matrices 
respectively. The three cases are labelled by the parameter /3, with (3 = 1 for orthogonal 
symmetry, /3 = 2 for unitary symmetry and (3 = 4 for unitary symplectic symmetry. One 
reason for this labelling is that the joint probability density function of 2 x 2 matrices with 
these symmetries is generically proportional to |Ai — A2|'^ for near degenerate levels. 

Regarding the region of the eigenvalue spectrum being scaled, depending on the details 
of the ensemble under consideration the eigenvalue density P(i){X) can exhibit various 
asymptotic forms. Our interest is when 

P(i)(A) = for A<0, ^«^^)~^;^ A^oo. (1.1) 

This corresponds to the so-called hard edge scaled state. As a concrete example of a ran- 
dom matrix ensemble which exhibits this behaviour, consider an nxp {n > p) real Gaussian 
matrix X, each entry independently distributed as a standard normal, and form the cor- 
responding covariance matrix X^X. The joint eigenvalue probability density function is 
then proportional to 

flxf/'e-^''/' n l^'^'-^.l'' ^'>0' (1-2) 

1=1 l<j<k<p 

with a = n—p+2/(3 — l and (3 = 1. The same construction, but with X a standard complex 
Gaussian matrix gives f 1 1.2 1) in the case /3 = 2. The probability density function fll.2p for 
general parameters will be referred to as ME^^Ar(A"''/^e~^'*'/^). For general (3 > 0, fixed 
a > 0, and with A; H- \i/4p, (/ = 1, . . . p — >■ oo, the eigenvalues in the neighbourhood 
of A = have spacing of order unity, and moreover exhibit the asymptotic behaviour (II. ip 

We therefore see that the hard edge scaled state is specified by the symmetry parameter 
/9, and the power law exponent a/3/2 (which represents the microscopic repulsion from the 
origin) in (11.21) . The notation E^^^'^{n\ (0, s); a/3/2) will be used to refer to the probability 
that in the hard edge state so specified, there are exactly n eigenvalues in the interval (0, s). 
Similarly, the notation £'^(n; (0, s); ME^ 7v(A°'^^^e''^'^/^)) will be used for the probability 
that in the ensemble specified by (11.21) . there are n eigenvalues in (0, s). 

The most prominent application in physics of the hard edge state is to the study of 
eigenvalues in lattice QCD. This comes about by the use of the matrix structure 



H 



pxp 



:i.3) 



to model the effective Hamiltonian for given topological anomaly u = n — p, and with 
elements of X independent standard Gaussians which are real (/3 = 1), complex (/3 = 2) 
or real quaternion (/3 = 4). The square of the non-zero eigenvalues of (II. 3p are equal 
to the eigenvalues of (11.20 . when the latter corresponds to the eigenvalue PDF of (11.20 . 
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Significantly, the exact distribution of tlie k-th smallest eigenvalue, p^^'^'^{k; {0, s); l3a/2) 
say, which is related to (0, s); (3a/2)} by 

, k 

p'^^^'^ik; (0, s); Pa/2) = Ef^\l- (0, s); /3a/2), 

1=0 

can be compared against data from lattice QCD simulations [T2] . 

Very recently, the explicit s ^ oc asymptotic form of E^'^^'^{n\ (0, s); (3a/2) for general 
/3 > 0, but requiring /3a/2 G Z>o, has been derived up to an error term which goes to 
zero [23]. The derivation relies on an unproved conjecture relating to the asymptotics 
of a certain generalized hypergeometric function. Such multivariable special functions, 
introduced into random matrix theory by Constantine and Muirhead in the case /3 = 1 in 
the 60's and 70's (see [30]), are finding their way into a number of recent works relating to 
asymptotics of eigenvalue distributions [3 [33 [2S1 [ISl El [IS] • The result of [23] also assumes 
knowledge of the corresponding asymptotics of £'^'^'''^(0; {0, s); (3a/2). The latter was first 
studied using generalized hypergeometric functions in [H], and the result 

1/9, /I \ a(/3a/2+l)/4-/3a/4 / / 1 \\ 

4-^^(0; (0, .); /3a/2) ^^^e-/^^/«+/^-^^V2Q ' ^^^/^{l + o{^)) , (1.4) 

where 

= 2(-^/^)'^(i-)"^'' n r(2,7/3), (1.5) 

was rigorously estabhshed for f3a/2 G Z>o, 2/13 E Z+. Soon after Dyson's log-Coulomb gas 
method was used to derive (but not prove) this same asymptotic form, up to the explicit 
form of the constant, for general paramaters [5]. In the case /3 = 2, (11. 5p has been proved 
for general Pa/2 > —1 in [6], and for |/3a/2| < 1 in [TT]. Very recently [33], stochastic 
differential equation methods based on tridiagonal matrices realizing the eigenvalue PDF 
(11. 2p have been used to prove (II. 4p for general /3 > 0, Pa/2 > —1, up to the explicit form 
of the constant. 

These developments motivate us to further consider the asymptotics of the hard edge 
gap probability E^'^'^'^^n; {0, s); Pa/2), both for n = and general n G Z"*". In Section 2 
we summarize previous work on this problem, as has been deduced using the method of 
generalized hypergeometric functions. In particular, the theory leading to (ll.4p is revised. 
The aim of Section 3 is to further develop this theory, so that the form (ll.4p can be 
rigorously established for Pa/2 E Z>o and general P > 0. In Section 4 we make use of a 
Gaussian fluctuation formula for linear statistics to determine a large deviation formula for 
Ei3{n; {0,x); ME^_Ar(A'''^/^e^^^/^)), first in the case n = 0, then for general n. Substituting 
X = s/ (4iV)^, then taking N ^ oo, these large deviation formulas are shown to scale to the 
asymptotic formulas for E^^^'^{n; {0, s); Pa/2) of Section 2, and furthermore give formulas 
which hold for general Pa/2 and P > 0. This latter point requires a suitable continuation of 
the product in (II. 5p beyond positive integer values of its upper terminal, which is achieved 
by making use of the Barnes double gamma function. As a check on these formulas. 
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an asymptotic functional equation known from [16], which relates the asymptotics of the 
spacing for a particular n, a) to those for parameters (4//3, /3(n+l)/2 — 1, /3(a — 2)/2 + 2), 
is shown to hold true. 



2 Generalized hypergeometric function expressions 

It turns out that for /3a/2 G Z+, the gap probability Ep{n; (0, s); ME^,Ar(A"'^/2e^^^/2)) can 
be recognized as an integral representation of certain hypergeometric functions based on 
Jack polynomials. Moreover, the hard edge scaling limit corresponds to a confluence limit 
of the hypergeometric function in question, giving back another hypergeometric function 
for which there is also a known integral representation, now as a /3a/2-dimensional integral 
and further requiring 2//3 G Z+ |14] . 

The class of generalized hypergeometric functions in question can be defined by the 
series 

pF^ \ai, . . . , Qp] 6i, . . . , bq] Xi, . . . , Xm) '■= ^ ] T/diTjTTW ^('^l' ' ' ' ' •^"^)' (^-l) 

K II' PiJk • • • [bqlK 

Here the sum is over all partitions > ^2 > ■ ■ ■ > > of non-negative integers, 
\k\ := Yl^=i generalized Pochhammer symbol [m]!"^ is defined by 



n(a--(j-l))^ , {a)k = a{a + l)---{a + k-l). (2.2) 

The function ci"'' (xi, . . . , Xm) is proportional to the Jack symmetric polynomial (see e.g. [T71 
§12.6]), and as such is a homogeneous symmetric polynomial of degree \k\. For m = 1, 
Ci"'* = x'^i and (12. ip reduces to the classical definition of pFq in one variable. Like pFq in 
one variable, pFq"^ exhibits the confluence property 

lim pF^")(ai, . . . , a^; 6i, . . . , 6^; — , . . . , — ) = p_iF^°^(ai, . . . ,ap_i; 6i, . . . ,bq; Xi, . . .,Xm)- 

Op^OO ^ V ttp ttp / ^ 

(2.3) 

In some cases there are explicit integral formulas for the generalized hypergeometric 
functions. One example is 

^(/3/2) 



,F^P''>{-N;a + 2m/^;tu...,tJ 



CfSa/2+m,l3 

/■oo poo ^ 

/ dx,--- dx^ll[xf'e-^'^/'llix,-t,)) H \x,-x/ (2.4) 

.7=1 fc=l l<j<k<N 



where C/3a/2+m,/3,m is equal to the integral with tk = 0, {k = l,...,m). This integral 
formula, obtained explicitly in [13], follows as a corollary of a similar integral formula for 
2Ff due to Kaneko [25]. 
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On the other hand the gap probabihty -E'/3(0; (0, s); ME/3_7v(A"''''^e~'^'^/^)) is obtained 
from (11.21) by integrating each A; over (0, s). Changing variables A/ H- A/ + s allows this 
gap probability to be recognized in terms of the integral in (12 ■4p . provided /3a/2 G Z>o, 
thus showing that in this circumstance [1] 



E0{O; (0, s); MEisA>^''^^'e-'''^')) = e-^''''\F[^"'\-N- a; {-sf (2.5) 

Here, in the argument of iF[^^'^\ the notation {uY means u repeated r times. Similarly, 
since by definition 



_ {N)n Cj3a/2,P,N 



\yj-yi\ 



n\ Ci3a/2,l3,N+n Jo 



l=\ l<j<l<n 
_ _ AT n 

X 

j=l l'=l i<j<k<N 



/3 



(2.6) 



we see from (I23D that for /3a/2 G Z>o, /3 G Z+ [23] 
Ep{n; (0,.);ME^,^+„(A'^'^/V^^/2)) 



n! 



C/3a/2,/3,7V+n l<j<l<n 

x,Fi^/'\-N- a + 2n; {-sT^/^ (-yO^ • • • , (-yn)"^)- (2.7) 

To proceed from (12. 5p and (12. 7p to the computation of £'^'^'''^(0; (0, s);/3a/2) requires 
that 

s ^ s/m, N ^oo. (2.8) 
Using the confluence (12. 3p it is a simple exercise to deduce from (12. 5 p that for /3a/2 G Z>o 



Ej'^'-'^(0; (0, s); /3a/2) = e"'^^/^ o^^i'^^'V; (s/4)'^'^/2)_ (2.9) 

Similarly, using (12. 3 p and (12. 7p . together with knowledge of the explicit gamma function 
functional form of Ca,i5,N [13 Prop. 4.7.3] one has that for /3a/2 G Z>o, /3 G Z+ [23] 

Ej^'"'^(n; (0, s)- I3al2) = Ap{n, a)s"+(/3/2)n(n+a-i)g-/3s/8 

/I /.I " 

X ,Fi^'^\- a + 2n; {s/Af^l^ {sy./Af, {syjif), (2.10) 
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where 



Next, we want to revise how (12. 9 P and also fl2.10p can be used for purposes of computing 
the s oo asymptotics. In relation to (12. 9 p this follows from a /9a/ 2- dimensional integral 
formula for iF[^^'^^ in (12. 9p . In addition to the requirement f3a/2 G Z>o, the integral 
formula also requires that 2//3 G 'L'^ . Under these assumptions, we then have [1] 

, /In (-l+2//3)/3a/2 / 1 n a/3/2 

4->;(0,.);/3a/2) = i?.,,e-^^/«Q^ 



/ TT gsV2eose,g^(-l+2//3)^^, TT ^^^e, _ ^^e,\^/P^Q^ . . . dOp,/,, (2.12) 



where 

_"^^ r(i + 2//3)r(2j7/3) 
^-^'^ - 11 r(i + 2j7/3) • ^2-^^) 

It was by using Laplace's asymptotic method to this integral that (II. 4p was derived. In 
[13] the asymptotic form (II. 4p was conjectured to hold true for general /3 > 0, and general 
/3a/2 > —1 (the latter subject to an appropriate interpretation of the product in (11.50 . 
identified in [23] as relating to the Barnes double gamma function; see (14.131) below). 

For the generalized hypergeometric function appearing in (12.101) there is no explicit 
integral form analogous to (I2.12p . Nonetheless a conjectured asymptotic form is available, 
which states that for s — j- oo and yi, ...,?/„ ~ 1 [231 eq- (3-45)], 

oFi^/'\- c- (s/4)^'^/2, {sy^/4y, {syjAY) 

= oF[^''\j,C- (s/4)/^('^+2n)/2)g/3V5E?.,(l-%)/2(^l + O(-ij)), (2.14) 

and this used in (I2.10p gives a conjecture for the s — )■ oo asymptotics of E^^'^'^iji] (0, s); f3a/2) 
Conj. 8]. 



Conjecture 1. For s — > oo we have 

E^-d(n;(0,s);/3a/2) ^^^^ f r ^ fn^ nax 



r^:/% (n) exp ( - /3 { - v/^n + + ^) log . V2 } ) (^1 + o 



i?hard(0; (0,s);/3a/2) 
where in the case [5n G Z>o 

2-{a+n)Pn . ^ . „(a+„-l)/3/2 ^ Y{a + 2j7/3) RJIo + (j + l)/3/2) 



(2.15) 



.hard / 
'/3a/2,/3l 



(2^)1/, n?:;'r(l + 0+a)/3/2)' 



j=l ^ ' llj=n 
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3 Integral formula for £;^^^^(0; (0, s); /3a/2) with pa/2 G 
Z+, /3 > 

We would like to generalize fl2.12p so that the restriction 2//5 G Z"*" can be removed. 
For this we return to the finite gap probability formula ( I2.5p . which has no such re- 
striction. Thus with j3a/2 G Z>o, but /3 > general, we can use a /3a/2-dimensional 
integral form of iF[^^'^\—N]a]{—sY'^/'^), deduced in turn from an integral formula for 
2F[^''^\-N,h]C] {sf"/^) [HI Exercises 13.1 Q4(i)] to deduce from ([23]) the following for- 
mula. 

Proposition 1. We have 

Ep{Q- (0, s)- MY.p^^iX'^P/^eP'''''^)) = e-^""''^ ^ 



M;3a/2(iV,-l + 2//3,2//3) 

1/2 1-1/2 



X 



J~\I1 J -111 ^ 



2-Kixi \ (- l+2//3+Ar) se^'f'^i 



where M„(a, ,6, c) denotes the Morris integral J77| , (Eq. 4-4)] > which is evaluated as a 
product of gamma functions. 

Proof. We indeed proceed as in pTf Exercises 13.1 Q4(i)]. Thus we begin with the 
formula [HI Eq. (13.11)] 



1 pl/2 pl/2 

— / dxi--- dxn ne™'('^-'')|l + e2"^^T+^(l + te^"^^')"' 

6, l/a) J_y2 J^l/2 fj[ 



X 

l<j<k<n 



Replacing t hj t/r and taking r — )■ oo gives 

•1/2 /.1/2 



\^2niXh _^2niXj \ 2/a 



— /" dxi--- [ TTe™'('^-^)|l+e2™t+'e"*''"''' TT 

Mn{aAl I a) J -1/2 J -1/2 t\ l</i<n 

= ,F("H-6;a + l + (n-l)/a; (t)"), 

where on the RHS use has been made of fl2.3l) . Recalling fl2.5p . we obtain a /3a/2- 
dimensional integral formula for £"^(0; (0, s); ME^^7v(A"'^''^e''^'^/^)) by setting n = /3a/2, 
a = P/2, t = —s, b = N and a = 2//5 — 1. Simple manipulation of |1 -)- e^™' | then gives 

(EU). □ 
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Essentially the same integral formula f l3.ip . but restricted to 2//3 G Z"^, was used in 
to compute the scaled limit (12.81) and thus derive (12.121) . The utility of the assumption 
2//3 G Z"*" is that the integrand can readily be rewritten to be analytic in 

Zk:=e^^''\ (fc = l,...,/3a/2) (3.3) 

except for a pole at the origin, thus allowing for the integration domain \zk\ = 1, (/c = 
1, . . . ,/3a/2), to be deformed. Deformation of the contours played a crucial role in the 
computing of the scaled limit (12.81) . We will now show that it is possible to deform the 
contours, and so compute (12.81) . without the need to assume 2//3 G Z"*". 

Proposition 2. Let C he the contour which starts at the origin, runs along the negative 
real axis in the bottom half plane to z = —1 — Oi, then along a counter clockwise circle to 
z = —1 + Oi, and finally back to the original long the negative real axis in the upper half 
plane. Let Ba^j^ be specified by h2.13\) . and assume /3a/2 G Z>o. We have 

, / 1 \ (-l+2//3)/3a/2 / 1 \ a/3/2 

^^(0; (0,.);/3a/2) = B^,,e~^^'^ {^-^' {£) 

(3.4) 

Proof. We begin by noting that the integrand in (13. ip is symmetric in the variables {xj}, 
so we are free to choose the ordering 

- 1/2 < Xi < X2 < ■ ■ ■ < Xpa/2 < 1/2, (3.5) 

provided the integral is multiplied by (/3a/2)!. But with this ordering 

( 2sin7r(xfc — Xj) 

l<j<k<(Sa/2 l<j<k<Pa/2 



X 



n \e'-^-e'--^\'/^= n (2sin7r(x,-x,))''''', (3.6) 



which is an analytic function of each of the variables (13. 3p in the corresponding complex 
planes cut along the negative real axis. The factors nf=i^ Q2TTtxi{-i+2/p)^-^_^^-2mxi-j(-i+2/i3+N) 
share this property, which thus becomes a property of the integrand. 

Let CN,t denote the contour which starts at 2; = — 1, runs along the negative real axis 
in the bottom half plane to 2; = —N/ \ft — Oi, then along a counter clockwise circle to 
z = —N/\/t + Oi, and finally back along the negative real axis in the upper half plane to 
z = — 1 + Oi. The analyticity properties of the integrand just discussed allow us to deform 
the original unit circle contours in the variables (13. 3p to the contours Cat^j, provided an 
ordering equivalent to (13. 5 p is adopted. Furthermore, in the variables (13. 3p . the integrand 
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and measure in (13.11) reads 

Pa/ 2 



-Q gsV2(,^.+l/^^.)/2^{-l+2//3) ^ 



1 N -l+2//3+Ar 



X TT f(-. --))'''^---^^. (3.7) 

V Zk J 2'KiZi 2'KiZRal2 

l<j<k<l3a/2 " ^ ^ '^"'^ 

This being symmetric in {zj} we can ignore the ordering constraint, provided we remove 
the earher introduced factor of (/3a/2)!. 

Now change variables Zj H- {N/y/s)zj. We see that in the scahng hmit ( 12. Sp . (13. 7p 
becomes proportional to (l)("~i+2//3)/3a/2 times the integrand in (13. 4p . thus accounting for 
all s dependent terms in the latter. The computation of the proportionality relies on also 
computing the asymptotics of the Morris integral in (13. ip . But this is the very same task 
as already carried out in [H], and the factors as presented in (13. 4p result. Finally, we note 
that with this change of variables and in the — oo limit each contour C^^t becomes the 
contour C, with the integral being well defined thereon. □ 

We remark that in the case 2//3 G Z"*", the integrand no longer has a branch cut along 
the negative real axis. The contributions to the integral of the portions of C running along 
the latter therefore cancel, and after parameterizing the remaining unit circle, we reclaim 

(Eia). 

The computation of the large s asymptotic form of (13. 4p is essentially the same as 
done in [13] for (I2.12p . The reason is that method of stationary phase tells us that the 
maximum contribution to the integrand then comes from the neigbourhood of Zj = 1. 
Thus the portion of C running along the negative real axis plays no role in this limit. The 
working of [H] is therefore justified without the restriction 2//3 G and so we obtain 
(II. 4p . proved now for /3a/2 G Z>o, and general (3 > 0. 



4 Gaussian fluctuation formulas 

The first aim of this section is to derive the large asymptotic form of 

E^{0- (0,S);ME^,;v(A'^^/^e-2/^^^)). 

We will see that by scaling the corresponding asymptotic formula according to 

s ^ s/{4Nf, TV ^ oo, (4.1) 

we can reclaim (II. 4p . We will demonstrate too that an analogous procedure can be applied 
to compute the asymptotics of the scaling limit for there being n eigenvalues at {?/j}j=i,...,n 
in (0, s) in the ensemble ME(3,N+niX'''^^^e-^'^^^). 
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4.1 Strategy 

Our method relies on reformulating the problem as one of computing the asymptotic form 
of a moment of the characteristic polynomial for ME/3^Ar(e~^''^^). In two recent works 
by the author [211 120], the spectral densities for the Gaussian, Laguerre and Jacobi /3- 
ensembles have been similarly formulated, and the corresponding asymptotics computed 
by recognizing that such averages can be interpreted as the characteristic function for the 
linear statistic V{x) = ^f^i log — A;|. The significance of this is that for ME^^Ar(e~^'^^'^), 
which is the scaled Laguerre ensemble with a = 0, eigenvalue density p(i),Ar(t) supported 
to leading order on (0, 1), it is a known theorem [3] that 

n{x + Xiy\ ~ gCMivWg(ca(.))V2^ (42) 

i'=l / MEa,jv(e-2/3iVA) 

where, with v{t) := log(x + t),x>0 

I^n{v) = j pi^i)^N{t)v{t)dt (4.3) 



1 2 n f 1 1 \ 

0'k = -J V \^- + -cosOj cos k9 d9, 



(4.4) 



and the notation ~ means that in the N oo limit the ratio of the LHS and RHS 

tends to unity. In words (14. 3 p says the characteristic function for the linear statistic exhibits 
Gaussian fluctuations, with an explicit mean (which is of order A^), and an explicit variance 
(which is of order unity). 

We will now show proceed to give the details of how this formalism applies to computing 
the sought large asymptotic forms. 

4.2 Large N form of EfsiO] (0, s); ME^^^vlA^^/^e-^/?^^)) 

Our first task is to show that E^(0; (0, 5); ME/s^Ni^'^^^^e-^'^^^)) can be expressed in terms 
of the average on the LHS of (14. 2p . 

Lemma 1. Lei A^(ME/3^jv(x^"/^e~''^/^)) denote the normalization required to make lil.^) a 
probability density function, and let s = s/4N. We have 

= e-W2(4Ar)^^ Ar(ME^,^(e-^"/^)) / ^^/A 

^ ' Ar(ME^,^(x/5-/2e-/3-/2)) + / • ^^-^^ 



ME^,jv(e-2/3'VA) 
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Proof. A simple change of variables A; i— A; + s in the definition shows 
E^{0- (0,s);ME^,;v(A^"/^e-^^/2)) 



Ar(ME^,Ar(a;/5«/2e-/3-/2)) 



'=1 / ME^,jv(e-/3^/2) 



We see from the definitions that with s = s/AN. The result (14.51) now follows by noting 
that the change of variables A/ t— )■ 4A^A; in the average on the RHS implies 

N \ IN 

X{{s + hf^i' j = {mf^ / His + Xif^i' 

□ 

We know from [22l eqns. (6.21) and (6.22)] and |2ll Eq. (4.1) with N = M] that for 
ME/3,w(e-2^^^) the density on [0, 1] is such 

2tP(i),^(t^) = - + - i) (^^' ^^'0 + °(^)- ^^-^^ 

With this explicit form, we want to compute (14. 3p . The following integral evaluation is 
required. 

Lemma 2. Let x > 0. We have 

log \ix + t\{l- tyl'^ dt = Vx(x + 1) - X - log (2{Vx^ + 1 - ^x)) - ^. (4.7) 

Proof. This can be deduced by an appropriate analytic continuation of [21], Eq. (3.2)]. 

□ 

Corollary 1. We have 



/ log |s + t| P(i),Ar(t) rft 
^0 



= 2iv( v/I(5TI - S - log {2{VpTi - VS)) + log (4-8) 

Proof. This follows by first noting 

log|s + t|p(i),Ar(t)rft = J log|iA/I + t|(2tp(i),Ar(t^))(it, 

where on the RHS 2tp(i) 7v(^) is given by the RHS of (14.361) . extended to be an even function 
by the addition of (s{t + 1) - ^(t)). Now (gS]) can be read off from fOTj) . □ 
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For the variance as specified by (14 ■4p . a straightforward modification of the working 
which gave Eq. (4.4)] shows 

cr' = -| log (5(5 + 1)) - 1 log (2((5 + 1) - S^^)) . (4.9) 

Substituting (KE\i in (H^ and (US]) in (lO]) . then substituting into (H^ we read off 
the large iV asymptotic form of the average on the RHS of (I4.5p . 

Proposition 3. We have 

N \ 

exp (iV/3a(Vs(S + l) - 5 - log (2((V52 + 1 _ y|) j _ 



4 V/3 2 ' ^ 



s 



X exp ( - ^ log(s(s + 1))^/' + ^ log ^ j . (4.10) 

Recalling (14. 5p . it remains to compute the large forms of the ratio of normalizations. 
Each is an example of a limiting form of the Selberg integral, and so is given by a product 
of gamma functions (see e.g. [HJ Prop. 4.7.3]), 

AA(ME^,^(e-/^-/^)) ^ . ./2^ W TT ^(1+^'/^/^) ^4 11 ^ 

Ar(ME;3,^(x/3''/2e-/^^/2)) ' ^ r (a/3/2 + 1 +j/3/2)" ^' ^ 

In the case that a G Z>o, simple manipulation of the product and use of Stirling's formula 
shows that the large iV form of (14. lip is 

a-l 

j^-aNfil2^j^^l^yfia{a-l)l^^^j^^yal2^aN{il2 "Q ^ ^.^^2). (4.12) 

j=0 

To compute the large form of (14. lip without the assumption a G Z>o, we follow the 
lead of the recent works [U |3T] and introduce the Barnes double gamma function V2{z\ 1, r). 
This function is related to the usual gamma function through the two functional equations 

1 ^ r-/-V2 r(z/r) 1 ^ 1 r(z) 

r2(z + l;l,r) 72^ r2(z;l,r)' r2(2 + r;l,r) V2^r2(z;l,r)' ^" ^ 

and is normalized by requiring that hm2_>.o 2;r2(2;; l,r) = 1. A result of Shintani [31] gives 
that r2(2;; 1, r) can be written as an infinite product of gamma fuctions, 

r.(.; 1. r) H2.r''e.p { (^ - |) logr + 

oo „/ N 2 
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where 7 denotes Euler's constant. Now, as shown in it is easy to check from the 
underlying recurrence that /^/2('^ + 1) := 11^=0 ^(1 + /3j/2) can be written in terms of the 
Barnes double gamma function according to 

Mn + 1) = (2^)n/2^-(n-^n(i-^))/2. r(n)r(l + n/r) ^ ^ ^ 2//^. (4.15) 

i 2[n, i, T) 

Further use of (I4.13P shows that this simplifies to give 

f,/,{n) = (2vr)("+^)/V-("-^W2r-"/^ ^^^^^l^^^ ^ ^^^^^^ 

But for a G Z>o we have that the product in (14.121) is equal to //3/2(o)//3/2(^)///3/2(^ + a)- 
This product is uniquely determined by its recurrence in N, and its value 1 at = 0, and 
thus the form implied by fl4.15p persists without the assumption a G Z>o. Moreoever, the 
asymptotic form of r2(2; 1, r) [32] then tells us that the large form of fi^/2{N)/ffs/2{N+a) 
continues naturally off its a G Z>o form, and so can be read off (14.121) . Thus the latter 
without the assumption a G Z>o reads 

where //3/2(a) is specified by KWh . Substituting KTOh and fliTTI) in (gS]) with c = l3a/2 
the sought large deviation formula results. 

Corollary 2. Let s = s/4N, and let f 13/2(0') be specified by lj^4-16\ ). With s fixed and 
a G Z>o we have 

e-^2A^'^' (iV/3/2) 1)^/^ (vriV/?) -'^Z (a) 

xexp|iV/3a(^^/I(JTl)-s + log(^VFTT+v^)) ^ T ~ ^) ^ ^ I } 

X exp I - ^ log(5(§ + l))^/2 ^ log ^ ^- 1. (4.18) 

An idea that goes back to Dyson [9] is to now scale the large deviation formula, which 
at the hard edge requires (12. 8p . to deduce the large s expansion of £'^^''^(0; (0, s); /3a/2). 
A rigorous justification of this procedure, which plays an essential role in a recent study 
of asymptotics of the soft edge gap probabihty [2], requires that a uniform error bound 
be provided with the large deviation formula [27]. Unfortunately the present method does 
not provide us with an error estimate. Nonetheless we find that the hard edge scaling of 
f|i?T8|) does reclaim fOj) . 

Proposition 4. Setting s = s/{4:N)'^, then taking N — )■ 00, the RHS of ( f^.igp is equal to 
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where 

Aa^p = (/3/2)-^'^('^-i)/^(7r/3)-'^/22-'^'^/2+V/3/2(a). (4.20) 
Moreover, in the case a/3/2 G Z"*" we have 

Aa,p = r^A, (4.21) 

where '^^ given by U.5\) . and thus ^.19^ reclaims the large s asymptotic form of 

^^^^■'^(O; (0, s);/3a/2) as given by (T^. 

Proof. Obtaining (I4.19P from fl4.18p is an elementary calculation. The equality f l4.2ip 
can be established by verifying that both sides satisfy the same defining recurrence, making 
use of on the LHS. □ 

4.3 Probability density function for there being n eigenvalues at 

{yj}j=i,...,n in (0, s) 

Let p{yi, . . . ,yn] (0, s); ME^ 7v+n(2^^''^^e~^^/^) denote the probability density function that 
in the ensemble ME^ 7v+n(2;'^''^^e~^^/^) there are n eigenvalues at {yj}j=i,...,n in (0, s). The 
aim of this subsection is to compute the large form of this probability density function, 
with 

fixed, and then the scaling limit 

So = so/(4iV)2, 5, =s,/(4iV)2 (j = l,...,n), N ^ 00 (4.23) 
of this hard edge form. 

Lemma 3. In terms of the notation lli4-22\ ) we have 
piyu . . . , (0, s); MEp,^^^{x^'^/'e-^^/')) 

^ A/" (ME^,jv(e-^^/^)) {N)n f,.r]m(a/2+n) -ms/2 fV 0a/2 

Ar(ME;3,^+„(a;«/5/2e-/3-/2)) n\ ^ ^ AJ-^' 

/ N n \ 

l<j<k<n \l=l j=l I ME^,jv(e-2/3iv^) 

Proof. We essentially follow the strategy of the proof of Lemma [H □ 

Comparing the average in fl4.24p with (12.41) we see that in the case /3a/2 G Z>o, fi G Z>o, 
we have 

MmpAe-^""")) ! ^. , . „ 



A/'(ME«^(a;('^+2"W2e-/'^/2)) \ J-J- 

= (4Ar)-^^W2+n)^^^(/3/2)^_^. ^ ^ (-4iVSo)'^"/', {-^NS^Y , . . . , {-ANs.,Y). 

(4.25) 
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(4.26) 



Furthermore, we recognise the average as the form given for large N by the RHS of (14. 2 p 
with c = 1 and 

v{t) = y log(t + So) + /3 J]log(t + Sj). 

Thus our immediate task is to compute the mean and variance. 
Lemma 4. With v{t) given by ^.20^ , and pi,(Ar)(t) as implied by ^4-30^ we have 

' (5a 



IJn[v) 



2iv(^) (v/5o(5o + l) - So - log (2(^/5^71 - 

n 

+2N0 J2 (\/%(% + 1) - >i - '°g (2(\/3K^ - y% 



2/3 4/12 ^ 



l + 5o 



So 



+ /3^1oj 



1 + Si 



Proof. We proceed as in the proof of Corollary [TJ 
Lemma 5. Let Uk := -(25^ + 1) + 2(s^ + Sfc)^^^- VFe /iave 



(4.27) 
□ 



iaiv)r = 2/3 



2. 



^) 1ok(1-z/2) 



Uj) — a 



-2 5^ log(l-z/,,z/,j}. 



(4.28) 



Proof. In the notation of (14. 4p the first task is to compute 



2 r/3a, cos6' \ cos6' \^ , ^ ,^ 



According to [211 Lemma 2], with Uk as specified above 



j=0 



where Wj = a/2 for j = and ii'j = 1 for j = 1, . . . ,n. Squaring this and performing the 
sum as required in (14. 4 p gives (I4.28p . □ 

According to the RHS of (14.21) we have that the large asymptotic form of the average 
in (14.240 is given by exp(/i7v(f ) + (a(w))^/2), with Hn{v) specified by (I4.27P and (ct(w))^ by 
(I4.28p . We are particularly interested in the double scaling limit obtained by performing 
(I4.23P in this asymptotic form when multiplied by appropriate prefactors. 
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Proposition 5. Let the large N asymptotic form of the average l[4-24^ , denoted by anno- 
tating the average by jj^, be scaled according to lli4-23\ ). We have 



/4^^7V/3(a/2+n) Ar(ME;3,jv(e "/^)) / ( ^/3a/2 TT.g. ^ 



# 



■5'=1 ' MEa,jv(e-2/3JV 



n n / I ^ — 

xexp(-4(|) logyi^+5:iogV37 + <.El°g(^^^^)} 



i=i i=i 

(-2/; ^ log( ^!^ )). (4. 



X exp 



where, with /^/2 specified by ^.161 ), 

Aa,n,i3 = (/3/2)-('^/^)("+2"-i)("+2")(7r/3)-(2"+")/2/^/2(a + 2n - 1). (4.30) 
According to fOSD and (Q, the scaled limit K23\f of the LHS of flCT]) is equal to 
oFr)(.;a + 2n;(|r^(|f,...,(^)/^). (4.31) 

Thus we obtain, as a conjecture, the corresponding large {sj}j=o,...,n asymptotic form. 

Conjecture 2. For /an^e values of the arguments {sj}j=o,...,nj i^^e generalized hypergeomet- 
ric function ([J]^^ has the asymptotic form given by the RHS of \4-2^ , up to terms which 
vanish as the arguments approach infinity. 

In the case /3 = 4 the asymptotic expansion of (14.311) can be derived from a matrix 
integral representation not available for general parameters [29] . The expression implied by 
the result of [29] is in precise agreement with the conjecture. Furthermore, this asymptotic 
form is also consistent with the conjecture (12.141) . 

In keeping with the origin of (14.311) . and with the scaled variables as specified by (I4.22p 
and (I4.23p . we can deduce from (14.241) that 

s>(syi, . . . , syn^ (0, s)- ym^^N^x^'^'^e-^^'^)) 

~ A,(n, a)."+(^/^)"("+'^-i)e-^^/« n(l - y.^/' \[ \y, - y,f 



But 



XoFr^.; a + 2n; (^)^«^ {^^^f, . . . , {^^^f). 



Ep{n- (0,.);ME;3,^+„(A'^^/2e^^^/2)) 

= / rfyi ■ • ■ / rfy„,p(si/i, . . . , syn, (0, s)- ME^,;v+„,(x^'^/'e-^^/2)) 
Jo Jo 
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and thus, recalling the consistency of Conjecture [2] with fl2.14p . we reclaim f l2.15p . but now 
with requirement that f3n G Z>o relaxed by writing 

|j r(a + 2j7/3) _ T2{a + 2/ (3; 1,2/ (3) 



(27r)V2 r2(2n + a + 2//3;l,2//3) 



^2n^ /r+2na/r-n/r+n /^/^(2n + a) 

(2^ fi/r{a) 



r :- 



2/13. (4.32) 



A significant check on fl2.15p with the substitution fl4.32p . supplemented by fll.4p and 
f l4.2ip . is to verify that it satisfies the asymptotic functional equation [16] 



Ef'\n- (0, s/~sp)- Pa/2) ~ Eljf{P{n + l)/2 - 1; (0, s/~s,/p)- a - 2 + 4//3), (4.33) 



where S4//3(/3/2)^ = sp. Note that on the RHS the value for the number of particles in the 
gap is/3(?7, + l)/2 — 1 and thus not necessarily an integer. In addition to fl4.32p we should 
also rewrite the remaining products in fl2.16p according to 



n?:o r(l + (j + l)P/2) _ f^^,{n + l)f^/2{n + a) 
n?:;'r(l + (j + a)/3/2) M2n + a) 

where fp/2{n) is given by f l4.16p . Combining with f l4.32p then gives 



(4.34) 



2-(a+n)2n/r ^n{a+n)/T+n /^/^(n + l)/i/^(n + O) 



{2txY fyria) 



2-{a+n)2n/r ^n(a+n)/r-l/2 ^ 2 - r)/i/^(r2 + a) 

r((n + l)/r) (27r)-+(i-)/2 J^^^) ' 



r = 2/(3, (4.35) 



where the second line follows from use of the definition fl4.16p and the functional properties 

dm. 

We showed in [21] that fl4.33p is satisfied up to the constant term in the corresponding 
asymptotic expansions of both sides. With §4/^ = 1, the constant term on the LHS of 
(^M is 

^mI^ r^^,";^,^(n)(/3/2)"(^"/2+l)/2-/3a/2+(n2+na)/3/2^ ^4^3g) 

while on the RHS the constant term is 

ra-2%M, r^-2%MM^ + l)/2 - !)• (4.37) 
Use of (14:351) and fOT]) shows that IK36\i is equal to 

2^ — {a+n)2n/T—a/T+a ^—1/2 

r((n + l)/r) ^2nr^a/2Hi~ry2 Mn + 1 - r)A/.(n + a). (4.38) 
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To verify the equality between f l4.37p and f l4.38p . we require the fact that the Barnes 
double gamma function has the inversion property 



T^in; 1, r) = ^-(l+"V2r)+n(l+r)/2rr^^!^. ^4^39) 

r r 



Recalling fl4.15p this implies 



(2Tr){a-l)/2 (l-a)/2 _ 1 

Ma) = y^-;^U— + 1). (4.40) 



Making use of this in f l4.37p shows the latter is equal to 

^-a' (a' -1)t /2-(a' /2+n')-n' [a' +n')t2-{a' +n')2n' t-a' t+a' 

r(n' + l)(27r)"'+"'/2 



/iK + l)/,(n' + a'), (4.41) 



where a' = [a — 2)t + 2, n' = [n + l)t — 1 and t = 1/r = /3/2. The expressions fl4.35p 
and fl4.2ip reveal that this is precisely the RHS of fl4.37p . thus verifying that the functional 
equation ( ]4.33p is a property of the asymptotic expansion of E^^^'^{n; (0, s); (3a/2) as implied 
by (O), (EISD, g35]) and I^M- 
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